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Abstract 

The long standing problem of inexplicably short spin relaxation in carbon nanotubes (CNTs) is 
examined. The curvature-mediated spin-orbital interaction is shown to induce fluctuating electron 
spin precession causing efficient relaxation in a manner analogous to the Dyakonov-Perel mecha- 
nism. Our calculation estimates longitudinal (spin-flip) and transversal (decoherence) relaxation 
times as short as 150 ps and 110 ps at room temperature, respectively, along with a pronounced 
anisotropic dependence. Interference of electrons originating from different valleys can lead to even 
faster dephasing. The results can help clarify the measured data, resolving discrepancies in the 
literature. 

PACS numbers: 85.35.Kt, 85.75.-d, 81.07.Vb, 75.76.+j 
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Spin-related effects in carbon nanotubes (CNTs) have attracted much attention in recent 
years Initially, it was assumed that a weak spin-orbital interaction (SOI) in CNTs, due to 
p^-electrons, would lead to a long spin relaxation (and spin coherence length) compared with 
typical semiconductors. However, numerous experimental measurements of spin transport 
and electron spin resonance in CNTs present very diparate results for the magnitude of 
the spin relaxation time r<j. The processing possibilities for spin transport impose spin 
conservation at least for the dwell time of the electrons in CNT providing a lower limit 
for Tg. Such approach, however, does not guarantee the direct derivation of r^, moreover, 
association of with time of electron drift between ferromagnetic contacts leads to enormous 
discrepancy with spin relaxation time measured by other methods. For example, a drift time 
around a ps attributes to magnetoresistance (MR) measurements in Refs. y and 14, while 
the electron spin resonance indicates much longer ^ 5 ns at temperature T = 300 K jsL 
Sufficient progress in the analysis of spin transport in CNT has been provided in Ref. |6 
where both = 60 ns and = 30 ns were extracted from MR measurements at T < 120 
K. 

A theoretical approach to spin relaxation based on hyperfine interactk)n failed to 



explain the short spin relaxation in CNT, even through a recent study [8| corrected the 
hyperfine interaction constant by two orders of magnitude, still making the spin relaxation 
longer than 100 jjs. While the carbon intra-atomic SOI has been evaluated as 12 meV j^, 
higher orders of perturbation theory result in extremely weak spin-orbital coupling (~ 1/ieV) 
for planar graphene 10|. Thus there is a negligible effect ori spin relaxation in CNTs. 
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Recently such an assumption was revised theoretically 10l-ll4| and experimentally 15 1 
because graphene curvature in CNTs produces mixing of n and a electron states that can 
drastically enhance SOI. Qualitatively, it can be viewed as electron spin in a magnetic field 
stemed from a clockwise or counterclockwise circular motion around the CNT's circumfer- 
ence. States with opposite electron motion originate from the two different valleys K and K' 



The finding of so strong SOI was shown to be responsible for an efficient mechanism of 



electron spin-lattice relaxation in CNT quantum dots at low temperatures 16[. The actual 
mechanism of itinerate electron spin relaxation in CNTs at room temperature still presents 
a very intriguing problem for carbon based spintronics. 

We begin with a few essential definitions concerning CNTs [17]. The CNT is conveniently 
imagined as a spiral graphite sheet (graphene) rolled along the chiral vector Ch (see Fig. 1). 
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Here Ch = naSL + n{,h, where a = a(l, 0) and b = a(l/2, \/3/2) are the graphene lattice unit 
vectors with a = 0.246 nm and Ua and rit are integers, which characterize the geometry of 
a particular CNT. The slope of Ch is defined by tan 9 = y/Srib/ {2na + rih) and the diameter 
of the CNT is given by d = \Ch \ /^r. The energy band structure of graphene possesses two 
non-equivalent valleys with Dirac-like dispersion law in the vicinity of Fermi level [l8|. They 
are located at the K = |^ (l, a/s) and K' = |^ "\/3) corner points of the first Brillouin 
zone, which will be labeled by A = 1, —1 respectively. 

The graphene two-valley band structure projects onto the CNT one so that in the vicinity 
of each valley the kP Hamiltonian takes the form Hx = hvF{\>im^i + ka2) where the 
Fermi velocity is fi;' = 8 ■ 10"^ cm/s and Pauli matrixes di and cr2 are defined over the 
sublattice electronic states A and B. The wave vector k with respect to (A = 1) or 
K' (A = —1) point is directed along the principal axis C, of CNT. The rolling along the 
perpendicular direction p (Fig. 1) imposes a quantization of the electron momentum in 
circumcircular direction fixing the wave numbers x^n = 2(m — Az//3)/(i with integer magnetic 
quantum number m and v (from the set 1, 0, -1) so that 2na + ni, + v becomes divisible by 3. 
If the CNT is subjected to a magnetic field Bq the Aharonov-Bohm magnetic flux passing 
through the CNT cross section (pAB = -Bo cos a ikP/A modifies the quantization condition 



x^^A = 2(m - Az//3 + </)ab/</'o)/c?, (1) 

where a is an angle between Bq and CNT axis, 0o = ^nficj |go| is the flux quantum, go the 
electron charge. 

The eigenvalues of Hamiltonian li\ describe the CNT electronic spectrum without spin 
structure, 



£:fc,r«,A = ±hVF^J xl^^x + (2) 

where -|- (-) corresponds to the conduction (valence) band. The CNTs with z/ = ±1 possess 
the semiconducting spectrum, with bandgap Eg = Ahvp/Sd even at zero magnetic field, with 
potential applications in spintronics We assume u = 1 in further numerical calculations. 

When a spin-orbital interaction is incorporated in the Hamiltonian Hx, the SOI mediated 
by CNT curvature which takes the simplest form in CNT co-ordinates ^, r], ( (Fig. 1) as: 

Hso = Ao\ao2Sc + Ai\ai2S^ - iA2a2{S+e''^ - ^.e"'^), (3) 
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where ao is 2 x 2 identical matrix, S± = ± iSr,, f = ^p/d, < p < \Ch\- In ([3]) the 
spin-orbital constants are proportional to CNT curvature, Aq = (So/d) cos3^^ , Ai^2 = Si,2/d 
where the parameter 6i = —0.19 meV-nm has been measured [3| and the ratios 6q/6i = 4.5 



and 62/ 6i = —1.4 can be estimated theoretically 10l-ll2|. The Aq depends on CNT chirality 
so that Aq is maximum in zigzag CNTs {6 = 0) and minimum (Aq = 0) in armchair ones 
{6 = 7T/Q) [l2|. 

These estimates show that the SOI may be treated as a perturbation for actual electronic 
energies, i.e., |Aj| ^ |£fc,m,A|- In the first order, the quantum averaging of the Hso on the 
eigenvectors \k,m, A) of Hamiltonian H\ gives rise to the reduced Hamiltonian of SOI H'gQ 
of electron spin energy in an effective field B'^^q directed along the principal axis ( of the 
CNT, H'gQ = gpsBgo^c, where gpsB'^Q = 2A(Ao + Aicti). We take into account that 
{k, m, X \ao\ k, m, A) = 1, {k, m, A |a2e^*'''| k, m. A) = and 

ai = {k,m,X\ai\k,m,X) = T—^==- (4) 

V ^ + 

Here =F in Eq- & corresponds to conduction (valence) band. 

An important issue of Eq. (jl]) is the dependence of WgQ^ai) on k i.e. electronic energy 
[eq. ([2])]. This effect can be understood as a result of diminishing of electron velocity v± in 
the p-direction with increasing fcy along the CNT axis under constant total electron velocity 
vp, i.e. v±/vf = k±/ |k| = k±/ ^Jkj_ + fcy that appears as a factor in Eq. (jlj) at k\\ = k and 

k± = Km,\- 

It is convenient to introduce the mean values (o"i) = 

T.m I ^1 f k,m,\)dk/'£^^ J f{ek,m,x)dk and B so = {B'so) = 2X{Ao + Ai {ai))/gpB, 
where f{ek,m,\) is a thermal population factor for electron with energy ek,m,x- Correspond- 
ingly, the SOI can be separated into a steady part gpsBsoSc, and a fluctuating one VtSc^ 
with 



= 2AiA(ai - (ai)). (5) 

Fluctuations in k lead to fluctuations of the effective spin-orbital fleld and serve as a mecha- 
nism of spin relaxation in CNTs. Our analysis centers on the efficiency of this mechanism for 
CNTs in the presence of an external magnetic field Bq. This supplements the Hamiltonian 
Hso with a Zeeman interaction Hz = (^/U^BoS. The joint action determines a spin energy 
gUB^effS in an effective field Be// = B50 + Bq. 
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In the reference frame t], (), the Zeeman Hamiltonian takes the form Hz = 
fkjJz{sma + cos a S(^), tvjjz = qiibBq {g = 2 and fiB is the Bohr magneton) so that 
the frequency of spin precession ux = gusBeff /h in a total effective field Bgjj is [20 1 



uJz sina)^ + [ujz cos a + 2A(Ao + Ai {oi))/K\ 



(6) 



Spin relaxation is conveniently described in other co-ordinates X, Y and Z, where Z is a 
spin quantization axis defined by the Bg//. Correspondingly, evolution of Sz and Sx, Sy 
correlates with longitudinal (spin-flip) and transversal (phase) spin relaxation. The angle 
C(x {0 < ax < tt) of Be// slope to CNT axis defines 



cosctA = [ItjjJz cos a -|- 2A(Ao -I- Ai {ai))]/'JTjjJx- 



(7) 



In a second quantization representation, the operator of fluctuating field takes the form 



^ = Z]fc,mb(fc,m)nfc^m, nfc,m = al ,^ak,m, where a{, „^ and ak^m are creation and annihilation 
operators and the vector b(/c,m) = 2AiA (cti — (cti)) (— sina^, 0, cosq;a) is expressed in co- 
ordinates X, Y, and Z. 

Assuming a short correlation time Tp for the electronic correlation function {nk^m{t)nk.m) 
in comparison with the spin relaxation time ts, one can use the Markovian equations of spin 
evolution in form of expanded Bloch equations in which the spin mean values (S) = Trp(t)S 
(p(t) is a density matrix j21|). Then the spin polarization deviation from thermal equilibrium 
becomes P = 2((S) - (S)^) , where 2 (S)^ =2 lim^^oo (S) = (0,0, - tanfkux/2kBT). These 
equations take the simplest form in the X, Y , Z coordinate system: 

dPx 



dt 



ux-T; 



(8) 



where the 3x3 matrix ux represents the cross product of T^x= (0, 0, ux), cDaPa = ^axPa- 
Thus, cDa includes only two non-zero matrix elements uyx = —^xy = <^a, and cDaSq = 0. In 
this case the matrix of relaxation parameters can be reduced to 





/ p(A) 
^XX 





^XZ 




fA = 





p{A) 





(9) 




d(A) 

\ ^xz 





4^1 J 





The low symmetry of the CNT affected by an arbitrary magnetic field produces multiple 
spin-relaxation parameters -R^'^v, n^u = X, Y, Z [four in the case of Eq. ([9])], which cannot 



be immediately associated with longitudinal or transversal spin relaxation rates. In the 
most general form they can be represented in terms of Fourier transformation of correlation 
functions {Q^(t)Q^)^ at frequency u = Ux. For spin precession in a fluctuating field they can 
be evaluated in an approximation of a momentum relaxation time 2M, l22| that assumes 
explicit knowledge of the dependence on k. Instead, we use the average value = (Tj^T^) 
to evaluate the spin-relaxation characteristics while the specific dependence of on k can 



be taken into account by introducing the numerical factor, which is about Qi ^ 1.5^3 



22|. 



Based on a recent study [23| we assign Tp = 1 ps and assume Qi = 2 for further numerical 
calculations. 

Applying these approximations we find 

i?S = cos' ax^^^^Q.Aal; 

Rfz = sin^«AT,|^^il^giAa?; (10) 

p{A) _ d(A) , d{A) . 
rtyy — rcj^^ -|- -fx^^, 

^(A) ^(A) 

Rxl = tan ax + ^ cot ax, 

where Acrf = {af) — ((Xi)^ critically depends on the temperature so that Ao"i — )■ at 
T — > 0. Results of calculations of the relaxation rate parameters [Eq. (11 Oj ) as a function 
of temperature, magnetic field strength and direction, and CNT diameter are presented 
in Fig. 2. It can be noticed that there is not a single Rj^J , which significantly dominates 
over all ranges of T, Bq, and d. Thus the full set of rI^J is required to describe spin 
relaxation in CNTs. All R^^^J are strongly depended on temperature that controls energy 
dispersion. Moreover, the relaxation rate decreases with an increase of CNT diameter due to 
SOI reduction with decreasing curvature. The non-monotonic dependence on the external 
magnetic field strength and its direction depends upon the interplay between B^o and Bq 
that leads to suppression of R^^l and Rx^z "when a ^ or Bq ^ 0. 

Examples of numerical solutions to Eqs. flHlQfTOj) are shown in Fig. 3. In particular, inter- 
ference of spin polarizations from the two non-equivalent valleys can lead to the oscillation 
amplitude beating shown in Fig. 3(a) where beats alternate with nodes in 90 ps. When the 
magnetic field is perpendicular to the CNT axis (a = 90°), a fast polarization damping is 
found for longitudinal (r^ ~ 150 ps) and transversal (r^ ^ 110 ps) spin polarizations as seen 
in Fig. 3(b). A greater variation of spin relaxation exists in CNTs with different chirality 



placed in an arbitrary directed magnetic field. 

The developed theory can be applied to MR measurements which are usually carried 
out on a CNT with weak contacts to source/drain ferromagnets. Electron undergo multiple 
reflections from the contacts before exiting. In such a case the output signal depends on 
both spin relaxation time Tg and electron dwell time [&]. The experimental setup also 
provides for electron injection into the CNT with spin polarization along magnetization M5 
of ferromagnetic source and spin detection with a ferromagnetic drain polarized in either 
the same direction M/) = M5 or the opposite one, M'^ = — M5. Spin dependent output 
from each valley A/a is proportional to difference mPA — hi'Pa = 2mPA (m = Mz)/ |M£)|, 
m' = — m) where the spin polarization is averaged over ( Pa = 'P\{t)e~*f'^'^dt) 
provided that is longer than electron drift from source to drain. Pa can be found from 
Eq. dH]) by multiplication with e~*/^'* and subsequent integration over t. Assuming that the 
rate of intervalley transitions is less than = l/r^ , the total spin-depended output can 
be expressed as A/a/2 = FsRd, where 

F, = ciny2(Tx + IRd-uJx)-'in. (11) 

A 

Here c is a dimensional constant and I is the identity 3x3 matrix. In the case of cua = 0, 
Fa = I/ts, Eq. (ITTll reduces to Fs = c/{t^^ + t^^) which is in agreement with previous 
calculations [6\. In more general situations, the following relation holds 

r. = - 1„„ "^hll. (12) 

Rd-^o dRd 

If uj\ ^ Rj^J = 1/tIu), the last equation leads to = ^x'^z'^/^x'^zz- Thus, only 
longitudinal spin relaxation (i.e. spin-flip) can control the magnetoresistance measurements 
in a CNT at strong magnetic field. In Fig. 3c the calculated dependence of on the chiral 
angle 6 of the CNT is displayed. For example, ^ 220 ps at T =300 K, for the case of 
Bq =0.5 T, d =2 nm, m || Bq a = 90° and 6 = 26°. 

However spin relaxation is suppressed at small values of a. The experimental setup 
treated in Ref. {g] displays a small deviation of the magnetization direction from the CNT 
axis, with a ~ 10°. For such a case and with T = 120 K, Bq =0.1 T and d = 10 nm Eq. ( 1T2|) 
predicts 69 ns which approaches the measured =30 ns. This results indicates the 
high efficiency of the spin relaxation mechanism even through small symmetry breaking in 
the design of realistic devices. 



In conclusion, we show that the fluctuations of curvature-induced spin-orbital interac- 
tion that are associated with electron random motion in CNTs is responsible for the short 
spin relaxation. If a magnetic field preserves the uniaxial symmetry of CNT, the spin-flip 
relaxation is suppressed. However a small asymmetry may lead to visible manifestations of 
this precession change. These issues may stimulate further experiments to determine the 
optimal conditions for room temperature CNT spintronic device applications. 

This work was supported in part by the US Army Research Office, NSF, and the FCRP 
Center on Functional Engineered Nano Architectonics (FENA). 
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FIG. 1: Reciprocal positions of CNT and coordinate systems. Insert: the lattice structure of 
graphene. 
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CNT diameter (nm) Tilting a of (degrees) 

FIG. 2: Rxx (solid lines), Rzz (dashed lines) and Rxz (doted lines) calculated as an average over 
valleys and the functions of (a) temperature [(na,n(,) = (15, 14), = 2T, a = '/r/2], (b) magnetic 
field strength [{na,nb) = (15, 14), T = 300 K, a = 7r/2 ], (c) the nanotube diameter [T = 300 K, 
5o = 2 T, a = 7r/2, 9 = 27°] and (d) the angle a between external magnetic field direction and 
CNT axis [(na, rib) = (17, 13), T = 300 K, = 2 T]. For convenience the right parts of the graphs 
are scaled in nanoseconds as a reciprocal values. 
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FIG. 3: (a) {na,nb) = (15,14), Bq = 0.2 T and a = 30° (b) (n^.n^) = (17,13), Bq = 0.05 T and 
a = 90°. Solid (dashed) lines corresponds to transversal (longitudinal P^) spin polarization. 
Thin horizontal line cuts the polarization curves at the level 1/e, which indicates the relaxation 
times 110 ps and 150 ps for and P(^. Insert (c) displays the dependence of relaxation time ts 
on chiral angle calculated with Eq. (jl2p under constant d = 10 nm, Bq = 2 T and a = 90°. 
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